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Abstract. We study a sudden interaction quench in the weak-coupling regime 
of the quantum sine-Gordon model. The real time dynamics of the bosonic mode 
occupation numbers is calculated using the flow equation method. While we 
cannot prove results for the asymptotic long time limit, we can establish the 
existence of an extended regime in time where the mode occupation numbers 
relax to twice their equilibrium values. This factor two indicates a non-equilibrium 
distribution and is a universal feature of weak interaction quenches. The weak- 
coupling quantum sine-Gordon model therefore turns out to be on the borderline 
between thcrmalization and non-thermalization. 



Sudden interaction quench in the quantum sine- Gordon model 



2 



1. Introduction 

The last years have witnessed an increasing interest in the dynamics of isolated 
quantum many-body systems. As it has happened to many other fields that were 
considered mainly academical in the past, the recent advances in experiments with 
ultracold gases confined in optical lattices have opened up the possibility of making 
real tests of the long time evolution of essentially isolated quantum systems. And 
following it, a pletora of unanswered questions have regained the attention of the 
scientific community. 

One important issue in this context is the question of thermalization in quantum 
many-body systems. Coupled to environments, the dynamics of quantum systems is 
known to lead to equilibrium states described by the thermal ensembles of Quantum 
Statistical Mechanics. But once the system is isolated and initialized in a highly non- 
thermal state, it remains unclear if the final, long-time state of the system can be 
described with one of those states, i.e., if the system has thermalized. Rigorously, 
it is easy to show that Quantum Mechanics does not allow evolution from a pure 
state to a thermal distribution, as unitary time evolution preserves the purity. 
However, to which extent thermal averages reproduce long-time quantum averages, 
and the conditions whereby this happens, can still be considered an open problem 

in [ii si in ei m m ei m nni nn H2i H3i m] - 

From this debate, new questions arose concerning the role that integrability plays 
in the long-time evolution of quantum systems. The dynamics of integrable systems is 
expected to be very constrained due to the large number of constants of motion, and 
hence ordinary thermalization should not be present in such systems. Experimental 
studies of systems near integrability point in this direction [15] . However, it has been 
argued that integrable systems could still relax to an ensemble described by quantum 
statistical mechanics if all the constants of motion are included as constraints [16 . 
In this respect the thermalization debate is still relevant even regarding integrable 
systems [T71IT8]. 

Similar to non-equilibrium classical statistical mechanics, the lack of a general 
framework to study non-equilibrium quantum problems makes it necessary to focus on 
specific models with the aim of extracting general features from them. However, long- 
time evolution of far-from-equilibrium systems is a very challenging topic itself, and 
only recently suitable techniques, both analytical and numerical have been developed. 
In this paper we employ the forward-backward scheme [21] 122] based on the flow 
equation method 23, 24 , which has turned out to be a realiable and powerful approach 
to solve Heisenberg equations of motion for operators. The main idea is to use the 
flow equation method to diagonalize the Hamiltonian in a controlled approximation, 
and then to study the time evolution problem in this diagonal basis. Hence real time 
evolution in this basis becomes simple and can be extended to long times without 
secular terms. All the difficulties of the problem are therefore encoded in the unitary 
transformation. 

In this paper we apply this approach to study the real time dynamics of the sine- 
Gordon model after a sudden quench of the interaction. Sudden interaction quenches 
are very interesting for studying non-equilibrium dynamics since they provide far-from- 
equilibrium initial states, but simplify the theoretical calculations since this initial 
state is simple. In optical lattices it is also possible to implement them experimentally 

nni- 

The quantum sine-Gordon model is a 1 + 1 dimensional scalar field theory with a 
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rich phase diagram. It is an integrable model, whose exact solution can be obtained by 
using the Bethe Ansatz ^U\. However, as usual this exact solution does not guarantee 
a simple calculation of the observables. Therefore it is interesting and necessary to 
implement approximate schemes like the forward-backward method. The choice of the 
sine-Gordon model is motivated by the fact that it is a paradigm for one dimensional 
translation-invariant interacting systems, with mappings that connect it to many other 
models. For example, it arises as the effective description of a system of interacting 
fermions in one dimension with backscattering, or for spin-1/2 quantum spin chains 
[25] . Sudden quenches have been already studied for one-dimensional fermions with 
density- density interactions, i.e., in the context of the Luttinger model |18j : this is 
an integrable model whose exact solution is quadratic once it is written in terms of 
bosonic excitations [35] . Likewise, the interaction quench to the Luther-Emery line of 
the sine-Gordon model leads to a quadratic model expressed in fermions |26j . However, 
in general such simple mappings are not possible for the sine-Gordon model: despite 
being integrable, it cannot be expressed as a quadratic Hamiltonian, and hence one 
can expect a redistribution of energy between different modes. 

The paper is organized as follows. In section 2 we introduce the main features 
of the model including its phase diagram and we describe its solution by the flow 
equation method. In section 3 we implement the forward-backward scheme to study 
the time evolution of the mode occupation operator in the weak-coupling region of 
the phase diagram. In section 4 the resulting expressions are employed to analyze the 
effect of a sudden quench of interactions: these are the central results of this work. 

2. The model and the flow equation solution 

The sine-Gordon model is an ubiquitous model widely studied in many different 
areas of physics. Its classical 1 + 1 dimensional version became very popular in 
the 1970s as it has non-perturbative solutions known as solitons [27]. Here we will 
be interested in its quantized counterpart, namely the quantum sine-Gordon model. 
As mentioned before the quantum sine-Gordon model is related to one-dimensional 
fermions with backscattering. There are many other similar mappings like to the 
one-dimensional Hubbard model near half-filling, the Coulomb gas problem, the two- 
dimensional classical X — Y model, and quantum spin chains [25] . Therefore the 
quantum sine-Gordon model is a natural and important setting for understanding 
quench dynamics. 

The Hamiltonian of the model is defined as follows: 

n = Jd x ^( x ) + l( d ^ + ^co smx ))) (i) 

where 4>{x) is a scalar bosonic field and H(x) its conjugate momentum field. In order 
to impose the quantum structure, they must satisfy the commutation relations: 

\n(x),<f>(y)]=-i8(x-y) (2) 

The Hamiltonian contains the parameter /3 and the coupling constant g, which define 
the phase diagram of the model. In the sequel we will usually use the parameter 
a 2 = /3 2 /47r, which will turn out to be directly related to the scaling dimension of the 
cos-interaction term. The rest of the parameters are used to regularize the theory: 
a is a lattice discretization parameter (its inverse l/a plays the role of an ultraviolet 
cutoff) and L is the system size (its inverse is the infrared cutoff). 
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Figure 1. Schematic phase diagram of the quantum sine-Gordon model 



A schematic picture of the phase diagram of the model is shown in Fig. Q] At 
a 2 = 1 (Thirring line) the Hamiltonian can be mapped to the non-interacting Thirring 
model [35] , whose relevant degrees of freedom are fermions [3D] > which can be identified 
with quantized solitons of the sine-Gordon equation. Away from the Thirring line 
these fermions experience an interaction: the region a 2 > 1 corresponds to repulsive 
quantum solitons, whereas in the region a 2 < 1 the interaction is attractive, which 
leads to bound states called breathers. 

Another interesting point of the phase diagram occurs near a 2 — 2. Here 
the system undergoes a Kosterlitz-Thouless continuous phase transition, which can 
be understood from the renormalization group equations for the flowing coupling 
constants 129 : 



dg 



dlog A 

da 2 
dlogA 



= (a 2 ~ 2)9 



4 2 

a g 



(3) 
(4) 



where A = 1/ v27ra is the ultraviolet cutoff. A graphical solution of these equations is 
shown in Fig. [5J For a 2 < 2 the coupling constant g flows to strong coupling, which 
signals the opening of a gap in the spectrum. This corresponds to the emergence of 
massive fermionic solitons as the appropriate low-energy degrees of freedom for the 
model. For a 2 > 2 we have the weak-coupling regime, where the coupling constant 
g flows to zero, and the relevant degrees of freedom are massless bosons. In this 
region, an approximate solution of the equations for fixed initial parameter a,® in the 
weak-coupling limit |<7o| "C 1 is simply: 



3(A) = 
a 2 (A) 



< A V 

A-0 



0{gl) 



(5) 
(6) 



In this paper we will be mainly interested in this so-defined weak-coupling limit 
where the flow of a 2 can be neglected to leading order. 
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Figure 2. Kosterlitz-Thouless like renormalization flow close to a 2 = 2 



Next we express the fields in terms of their bosonic modes: 



V4^ fc>0 Vk 



H(x) = -1= ^ e-^Vk (e- ikx {a{ k a r ,_ fe ) + e ife > ; , fe - <_ fc )\g) 

The various creation and annihilation operators obey the usual bosonic commutation 
relations: 



[ai,k,a\ k ,] = [a r - k ,al_ k ,] = 5 kk , — Vfc,fc' > 
' Ztt 

We use the Mandelstam vertex operators [3D] to construct the Hamiltonian: 

T7-/ \ „ / \ * —ah/2 1 / — ikx t ikx \\ 

Vl(a;x) =: exp(q } y e 7 y^( e a \.k~ e a Uk)) ■ 



(9) 



fe>0 



= (JL r 2 /2 exp(a £ e - Qfe /2 _L (e-tt-aj fc - e ^ ai ,0) 
27ra fc^S ^ 

K(«5 1) =: exp(a £ e- afe / 2 ^( e -^a r ,_ fc - e^a\ _ k )) : 
fc>0 vfc 



(10) 



(11) 



fc>0 



where : O := — (0|O|0) means normal ordering of the operator O with respect to 
the non-interacting ground state |31j . The Hamiltonian now reads: 

U = £ k{a\ k a Lk + al_ k a r ,-k) 

k>0 

2na 2{ L 

This Hamiltonian can be studied in its entire phase diagram using the flow equation 
approach [M as shown in Refs. [321 [S3]- The key idea behind this method is to 



+7±^(—) a I dx{Vi(a;x)V r (-a;x) + V r (a;x)Vi(-a;x)) (12) 
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successively apply infinitesimal unitary transformations that eventually diagonalize 
the Hamiltonian. The stable choice of such a transformation sequence requires 
energy scale separation similar to conventional renormalization approaches: as the 
transformation progresses, more and more energy-diagonal interaction matrix elements 
are eliminated. Such a scheme was proposed by Wegner [23] and independently 
by Glazek and Wilson [34] [35] . Wegner showed that a suitable infinitesimal 
transformation is obtained with the following canonical generator: 

r)(B) = [H (B),n int (B)} (13) 

where Hq is the diagonal and Hint is the interaction part of the Hamiltonian. B 
is the flow parameter that parametrizes the diagonalizing flow. It can be related to 
an energy scale A# — 1/vfl in analogy to the conventional RG scheme. The key 
difference is that the flow equation Kb corresponds to an energy difference that is 
being eliminated, whereas in a conventional RG scheme A correponds to the UV- 
cutoff, which is an absolute energy scale. 

The flow of the Hamiltonian is given by the following differential equation: 

^l = [r,(B),H(B)} (14) 

The methodological challenges come from the implementation of such a transformation 
and the integration of the ensuing differential equations. In most cases approximations 
are required in order to get a closed sets of equations. However, these approximations 
do not necessarily match those employed in other methods like perturbation theory, 
which, e.g., opens the possibility to access non-perturbative regimes using flow 
equations. For a comprehensive review of the flow equation method and its 
applications we refer the reader to Ref. [24] . 

Before studying the non-equilibrium dynamics of the sine-Gordon model, we first 
briefly review the flow equation solution of the equilibrium sine-Gordon model in order 
to make this paper self-contained. More details of this calculation can be found in 
Refs. [31 [33]. 

It turns out to be more convenient to work with Fourier transformed vertex 
operators: 

Vl (- a -k) = ^- f dxe- ikx Vi{-a]x) 
Air J 

Vi(a;k) =V?{-a;k) = / dxe ikx Vi(a;x) (15) 
V r (—a; k) = — !- / dxe- ikx V r (-a;x) 

V r (a;k) = Vj(-a;k) = ^- [ dxe ikx V r (a;x) (16) 

27T J 

Some relevant properties of these operators are summarized in Appendix A. In this 
representation the generator of the unitary transformation consists of two parts: 

rj(B) =?? (1) (B) +V {2) {B) 

(B) = 87r 2 ^^pu(p; B) (Vi(a; p)V r {—a; p) — h.c.) (17) 

p 

V W(B) = -V(£)$>W,-fc " <*r,-kai, k ) (18) 

k>0 
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u{p\B) = /0 s 2 (-^) * 



(27ra) 2V L 



32,32B^_ o3(B y 



« 2 (B) 



1 



(19) 



(20) 



47rr(a 2 - 1) 
(a 2 -2) 3 + 0( 3 2 ) 



(21) 
(22) 



4r(a 2 (£) 

The actual values for these coefficients are obtained by solving the flow equations for 
g{B) and a{B) [32l|33]: 

da 2 _ a 4 (g 2 + 0(g 3 )) 
dl 
dg 
dl 

with I = — | log(32£?/a 2 ). This flow succesfully describes the different regions of 
the quantum sine-Gordon phase diagram, from the weak-coupling regime close to the 
Kosterlitz-Thouless transition, where excitations are massless bosons, to the Thirring 
line at a 2 = 1 with massive solitonic excitations. 

Since we are interested in the weak-coupling limit (go small for fixed initial a) in 
this paper, we can neglect the flow of a(B) to leading order and identify a(B) with 
its initial value a. The effective diagonal Hamiltonian generated in the limit B = oo 
then has the following form: 

H(B = oo) = H + H diag (B = oo) 

Hq = H a l,k a l,k + a l-k a r-k) 

k>0 



(23) 
(24) 



H d iag(B = OO) = y^6Jfc(5 = OO) 



(25) 



fe>0 



x (P z (a; -fc)P/(a; -k) + P/(a; fc)P(a; fc) 
+P r f (Q;; -k)P r (a; -k) + P r (a; fc)P r t (a; fc)) 



Here 



pt(a;fc) 



2?r 
2tt 



r(a 2 ) 
r(a 2 ) 



L\k\ 
2tt 

2tt 



l-a 



2-, 1/2 



1/2 



(26) 

(27) 
oo) is 



are conveniently normalized Fourier transformed vertex operators. u>k(B 
given bv |52ll53]: 

u, k (B = oo) = -g 2 ft \ak\ 2 ^ (28) 



2r 2 (a 2 ) 



3. Implementation of the forward-backward scheme 

The flow equation approach to the quantum sine-Gordon model can be used to study 
the real time dynamics after a sudden interaction quench. The general idea proposed 
in [3T] consists of three steps. First, the observable is transformed into the diagonal 
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basis of the Hamiltonian, that is one carries out the same sequence of infinitesimal 
unitary transformations: 

=[„(£), 0(B)] (29) 

This is the so-called forward transformation of the observable, and in most cases it 
implies a very complicated structure of 0(B — oo). This observation is familiar from 
exact Bethe ansatz solutions. 

The advantage of working in the diagonal basis comes from the actual time 
evolution, which is much easier for a diagonal Hamiltonian. Later we will see that in 
our model the term Tiding from Eq. (|25|) can be neglected compared to Ho for not too 
long times. Hence, time evolution translates into phase factors and 0(B = oo,t) is 
obtained easily. 

The final step is the backward transformation, where the flow of the transformed 
and time-evolved operator back to the original basis is carried out. The result 0(t) is 
then an approximate solution to the Heisenberg equations of motion for the operator 
and it is straightforward to work out its expectation value with respect to the initial 
(non-equilibrium) state. 

This sequence of transformations constitutes the flow equation forward-backward 
scheme. It has been already used to study the real time dynamics of a Fermi liquid 
after a sudden interaction quench [51 [7] and the real time dynamics after a sudden 
quench in the ferromagnetic Kondo Model [H [9]. One of the main advantages of 
this approach is that it avoids the infamous problem of secular terms in perturbation 
theory: in time-dependent perturbative expansions, these can restrict the perturbative 
solution to time scales shorter than [coupling constant}^ 1 . In this sense, the forward- 
backward scheme is the quantum version of unitary perturbation theory in classical 
mechanics. 

In this paper we are mainly interested in the time evolution of the bosonic number 
operator after a sudden interaction quench. Hence, our first goal is the implementation 
of the forward-backward scheme for the creation/anhinilation operators. 

3.1. Forward transformation 

The forward transformation requires the solution of the flow equation 

^M = [m>ai , k (B)} (30) 

with the initial condition ai^{B = 0) = a, ( fe. Here i = l,r corresponds to left and 
right movers. In order to get closed equations we make an ansatz: 

a Ilfc (oo) = hf],(B)ai, k + h%{B)a\_ k (31) 

+4ir 2 ]T 4°(p; B) {Vi(a,p- k)V r {-a;p) - V r (a; -p)V5(-a; k - p)) 
p 

and likewise for the right movers. The ansatz is parametrized by various functions that 
must be calculated by working out the commutators in equation Q30p. In the weak- 
coupling phase this ansatz becomes exact in the infrared (low-energy) limit since the 
coupling constant g(B) flows to zero. 
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It is convenient to decompose the transformation in two stages given by T}W{B) 
and rj( 2 \B). The lowest order contribution from the first part of the generator yields 
the following flow equation: 



iUJ - -j= [h\l' k (B) + h%(B) ) pu(p; B) (32) 

Now we focus on the weak-coupling limit. Since the coupling constant g(B) flows to 
zero, we can use it as a perturbative parameter in the flow equations. It can be shown 
that in order to preserve the bosonic commutation relations during the flow 

[a itk (B),a\ k ,(B)}=6 id 5 ktk ,-^- (33) 

it is consistent to assume an expansion of the form: hf k (B) = 1 + 0(g 2 ) and 

h { ^ k = + 0(g 2 ) (see Appendices B and C). 

In the weak-coupling limit the flow of the coupling constants can be approximated 
by equations ([5]) and ©, which simplifies the integration of the differential equations. 
The result of carrying out the whole flow is: 

fij? (p; co) = P a2 - 3 ^ r(2 - y , (2paf) (34) 
where F a — ,„ &z^2 (^r L ) a ' ■ This result is valid for a 2 < 4, which is the region 

(2a) Q 

we are mainly interested in. 

Now let us discuss the effect of the second part of the generator r)( 2 \B). Due 
to the structure of this generator a different approach is possible. The complete 
infinitesimal transformation can be rewritten as: 

a^ k {B + dB) ~ a l>k (B) + [77(B), a^ k {B)]dB 

= e"' 2 ' (o,, fc (fl) + [rP{B),<H, k [B)]dB) ^ (35) 

The term in brackets corresponds to the transformation already worked out above. 
The advantage of this expression arises from the fact that we already know the effect 
of the exponentiated rf 2 '{B) on the bosons and vertex operators: 

e" (2)(B) a; , fce -" (2)(s) = a hk cosh(^(.B)) + a\ _ k sinh(V>(£)) (36) 

e' (i) < B »aJ ri e-' W ( B ' = a\ _ k cosh^B)) + a t , k sinh(^(B)) (37) 

e" (2)(B) ^(a;p)e^ <2)(iJ) ~Vi(a(l + 4>(B)ip) (38) 

e ri {2) {B)y^_ a . p y~r l ^{B) „ V r (-a(l + 1>{B)',p) (39) 

Hence the effect of the second part of the transformation is: 

e^^ aiik (B)e~^^ = 

(hfl(B) cosh(V>(B)) + h l rM {B) sinh(^(B))) a hk 
+ (hfl(B) sinh(V(£0) + hf k {B) cosh(^(B))) a\_ k 

p 

x (Vt(a (I + ^{B));p- k)V r (-a(l + tp(B)); p) 

- V r (a (1 + i>{B)) ; -p)Vi(-a(l + t/j(B)); k - p)) (40) 



Sudden interaction quench in the quantum sine- Gordon model 



10 



Carrying out the whole transformation in second order of the coupling constant yields: 



+4tt 2 4°(p: to ) - *)V r (-a;p) - V r (a; -p)VS(-a; fc - p)) 

p 

(41) 

Here we have used a decomposition derived in Appendix B: hf k (oo) ~ 1 — hf\.(oo) ~ 

1 — Sq^j,4 • We will later see that in the present order of the calculation the 
contribution coming from the second part of the transformation can be neglected 
since ^i(oo) oc g%. This is consistent with our assumption that we neglect the flow of 
a 2 (B) in the weak-coupling limit: neglecting the flow of a 2 (B) in fact just corresponds 
to neglecting the generator part rf 2 ^ . 

3.2. Time evolution in the diagonal basis 

The second step in the forward-backward scheme is the time evolution of the observable 
in the diagonal basis. Here, however, an additional approximation is required in order 
to solve the time evolution problem: in the diagonal Hamiltonian (|23|) only the bosonic 
kinetic term Hq is taken into account. We will later see that this approximation 
implies a maximum time scale up to which our calculation can be trusted. The time 
evolution dictated by is straightforward due to the simple transformation of the 
vertex operators: 

e m « t V l (-a,p)e- m « t = e^V^-a^) (42) 
e mot V r (-a,p)e- inot = e tpt V r (-a,p) (43) 
Therefore the time evolved anhinilation operator in the diagonal basis reads: 

ai,fc(oo,t) = ^1 - <?o 2 E4°J e "Vfc + (so 2 E4° ^ a l,-k 

+ 4Tr 2 e- ikt J2^\ P ;oo)e 2i P t (V l (a,p-k)V r (-a;p) 
p 

-V r (a;-p)Vi{-a;k-p)) (44) 

3.3. Backward transformation 

The final step of the transformation requires to undo the flow equation transformation 
for the time-evolved operator (|44|. This is straightforward due to the perturbative 
nature of the transformation. We make a general ansatz for the operator: 

a Lk (B,t) = hfl{t;B)e- lkt a Uk + h^(t; B)e ikt a\_ k 

+4^e- ikt J2n k l \p,t;B)e 2i v t (V l (a-,p-k)V r (-a;p) 
p 

-V r {a;-p)Vi{-a;k-p)) (45) 

with initial conditions hfj.(t; oo) = 1 - goJ2 p z p\ ^ k (t;oo) = g% J2 P z p ^ + VK 00 ) 
and Cl k l ^(p, t;oo) = Qj^(p;oo). The flow equations resemble those for the forward 
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transformation, since the contributions from the other functions in the ansatz follow 
directly via the bosonic commutation relation, we will only explicitly write down the 
flow equation for the function 0^ (p, t; B) : 

with the solution: 

^ ) (p;t;0) = ^ ) (p;^)(l- e - 2 ^) (47) 

By applying the first part of the generator the final time-evolved operator in second 
order of the rcnormalized coupling constant then reads: 

+4TT 2 e- M J2n k l) (p;oo)(e 2l P t - l)(Vi(a,p - k)V r (-a;p) 
p 

-V r {a;-p)Vi(-a;k-p)) (48) 
This is the main technical result of our paper, which can be used as a building block 
to study the time evolution of all other observables. 

3.4- Consistency check: Ground state energy in perturbation theory 

As a consistency check for the flow equation calculation we now compare the flow 

(2) 

equation result for the ground state energy E^ with second order perturbation theory. 
In fact it is sufficient to evalute the kinetic energy E K ' in the ground state since one 
can easily prove E^ Q — —E^\ Our goal is therefore to calculate 

E K ,o = (0\n K \0) =J2 k ((°Kfcl°) + (6K,fc|6» (49) 

k 

where |0) is the ground state of the interacting model. 

Within the flow equation formalism this is most conveniently evaluated in the 
diagonal basis with the forward transformed operators 

E Kfi = E k (<°Kfc(~Hfc(°°)l°> + <0|< fe (ooX, fe (oo)|0>) (50) 
k 

where we have used (Oln^/dO) = (0\ni^{B = oo)|0). Here |0) is the bosonic vacuum 
since this is trivially the ground state in the diagonal basis. If this calculation is carried 
out using our previous results from the flow equation formalism, the result actually 
does not coincide with perturbation theory. However, this is simply due to the fact 
that the flow equation calculation is a rcnormalized expansion, whereas conventional 
perturbation does not contain renormalization effects. In order to compare with 
perturbation theoryj we therefore artifically set g{B) = g n and find: 

feq _g 2 L T(2a 2 2) 
" a 2na T 2 (a 2 ) 

| And only for that reason since renormalization effects due to the running coupling constant are 
essential in the sine-Gordon model. 
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The same result can be obtained by working out the kinetic energy in second order 
perturbation theory: 

E { k\= dXy2\( n \^i\°)\ 2e ~ X ^ kHni ' h+nr ' k) ( 52 ) 

where \n) = (nj^j, ft.;,fc 2 , n r j^, n ri k 2 j ••)• By using the matrix elements of vextex 
operators given in Appendix A we arrive at the following expression: 

= 2L jSr I ^ jf ^ E **~ (ta)e ~" (53) 



After working out the integrals it can be shown that the flow equation result ([51]) is 
reproduced. 

4. Real time dynamics after a sudden interaction quench 

We now want to use the previous results to investigate the non-equilibrium dynamics 
of the quantum sine-Gordon model in the weak-coupling limit. Specifically, we study 
the mode occupation numbers after a sudden interaction quench: 

H{t) =H + Q(t)Hi (54) 

Since the system is prepared in the non- interacting ground state |0) (bosonic vacuum) 
of Hq for t < 0, the time-dependent occuptation number for left movers is simply 

(n, lfc )(t) = (0(t)|m, fc |0(t)) = (0|oJ fc (t)o,, fc (t)|0> (55) 

where we insert the time-evolved operators in the Heiseberg picture (J3HJ) • It is this 
conceptual simplicity which makes sudden interaction quenches very appealing for 
studying non-equilibrium problems. Notice that an identical result to (jH} can be 
obtained for right movers, and therefore we restrict ourselves to explicit expressions 
for left movers only. 

An expression for ([55)1 can be worked out readily from the result (|4"B")) obtained 
within the forward-backward scheme: 

a 2 g 2 L 4" 2 - 2 jn t 
'27ra(r(a 2 )) 2 
with the integral 

f r 00 f n 2 2 2 

I(ka,-)= / da;sin 2 (a;-)(r(2 , 4a; 2 ))V Q ~ 7 (fca - x) a - l e~ x 

a ha a 2 

In order to get this result, we have made use of the properties of vertex operators 

summarized in Appendix A and of the expression for f2? (p; oo) given in (|34[) . However, 

it is more convenient to express this result in terms of the equilibrium occupation 

numbers. Fortunately, the equilibrium occupation numbers follow directly from the 

flow equation calculation: 

Meq = (0Kfe|0) = (0\a[ k (B = oo)4 ik (B = oc)|0) 

n 2 2 r 4« 2 -2 

k 2na (1 (a' 2 )) z 

Here |0) denotes the interacting ground state, and: 

1 r°° m 2 2 2 

I EQ (ka) = - / dx(T(2 - ^-,4x 2 )) 2 x 3a - 7 (ka - x) a (58) 

4 Jka 2 



Mt)) - -^ — TW-^Hka, -) + O( 5o 4 ) (56) 
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Figure 3. Real time dynamics of the occupation number normalized to the 
equilibrium value for different values of a 2 after an interaction quench. All the 
curves correspond to the same momentum ka = 0.1. After a few oscillations the 
occupation numbers converge to twice their equilibrium values, see text. 



This results in the following compact expression for the ratio of non-equilibrium 
to equilibrium occupation number: 

{ni,k)EQ lEQ[ka) 

and likewise for right movers. The integrals in this expression must be computed 
numerically. Plots of these ratios are shown in Fig. [3] for different values of a 2 and 
fixed momentum, and in Fig. U for fixed a 2 and different momenta k a. 

One observes that the key phenomena after the quench are damped oscillations 
of the mode occupation on a time scale set by the lattice cutoff a. The asymptotic 
value of the mode occupation universally converges to twice its equilibrium value: 

(n I)fc (t -> oo)) = (2 + 0{gl)) {n Uk ) EQ (60) 

This can be understood easily by noticing that I^Q{ka) and I(ka,t/a) only differ by 
replacing a factor 1/4 by sin 2 (x^/a) in the integrand. Clearly the limit (nj fc(i — > oo)) 
just amounts to taking the time average over sin 2 (xt/a) in the integrand, which gives 
1/2 and therefore I(ka,t/a —> oo) = 2lEq(ka). 

Notice that (IBTfl) implies a non-thermal mode distribution function for the 
asymptotic state of this closed quantum system: the equilibrium system with nonzero 
temperature cannot reproduce this expression. 

5. Conclusions 



We have demonstrated that an interaction quench in the weak-coupling phase of the 
sine-Gordon model leads to an interesting dynamics that is quite different from a 




quench of the forward scattering only. Quenching the forward scattering does not 
induce any dynamics for the bosonic occupation numbers, which are in fact constants 
of motion in this case [18]. On the other hand, we have found that a quench of the 
backscattering term in the weak-coupling phase leads to a real time dynamics where 
the excitation energy of the quench is converted into bosonic mode occupations that 
oscillate on a time scale set by the ultraviolet cutoff and eventually reach twice their 
equilibrium values (I60|) . This factor 2 is universal for weak interaction quenches and 
has previously also been seen for the momentum distribution function in the non- 
equilibrium Hubbard model [6j [7] and for the magnetization of the non-equilibrium 
ferromagnetic Kondo model |9] . It occurs for quenches in quantum systems where 

i) second order perturbation theory is valid (at least up to a certain time scale) and 

ii) for observables like the mode occupation number operator which commute with Hq 
(for a proof and more details on the conditions see Ref. [7]). 

Since such general results are important for developing a better understanding 
of non-equilibrium quantum many-body systems in general, we need to critically re- 
examine the approximations in our calculation. Due to the weak-coupling behavior of 
the running coupling constant, the second order calculation presented here becomes 
more and more reliable in the infrared limit. Therefore higher order corrections to the 
universal factor 2 will vanish in the low-energy limit. However, we had to make the 
additional approximation to neglect the time evolution generated by Hdiagipo) from 
(|25[) . Therefore the result (|60|) can only be trusted up to the time scale 

r fc oc oj^{B = oo)(x g„ 2 k 3 ~ 2a2 (61) 
Whether our central result 



(n it k(t oo)) = 2(n iyk ) E Q 



(62) 
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really holds beyond the time scale or only in a time window a < f < cannot 
be answered based on our calculation. The interaction quench in the weak-coupling 
phase of the sine-Gordon model therefore occupies an interesting place between the 
ferromagnetic Kondo model, where the factor 2 is in fact asymptotically exact [HI IS], 
and the non-equilibrium Hubbard model in d > 2 dimensions, where the factor 2 
describes the prethermalization regime [6j[7] before the system eventually thermalizes. 
The integrability of the sine-Gordon model could make one expect that the non- 
equilibrium distribution function (|62|) remains stable for all times and does not 
approach a thermal limit form. However, the conserved quantities in the sine-Gordon 
model do not impose any obvious constraints on the dynamics of the momentum 
distribution function, unlike in the case of quenching the forward scattering 18J. 

Further studies of this question would be very worthwile, either numerical or 
based on an exact solution. We have seen that the weak-coupling quench in the sine- 
Gordon model is on the borderline between thermalization and non-thermalization 
seen through the eyes of the mode distribution function. This is similar to the role 
played by the celebrated Fermi-Pasta-Ulam problem for classical many-body systems 
[36] [37] . A better understanding of the weak-coupling quench in the quantum sine- 
Gordon model could be an important step in elucidating the fundamental question of 
thermalization in the quantum world. 
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Appendix A. Properties of vertex operators 

We summarize some important properties of vertex operators (k, k' > 0). More details 
can be found in Ref. 1331. 





V t (a; 



k)\0) = Vi{-a;k)\Q) = V r (a;k)\0) = V r {-a;-k)\Q) = 



(A.l) 



(o|vk 



a; k)Vi(a;k')\0) = (0\V r (a; k)V r (-a; k')\Q) 




(A.2) 



(0\Vi{a;k)Vi{-a;k')\Q) = (0\V r (-a;k)V r (a;k')\Q) 




2 



(A.3) 



(0\Vi(a;k)Vi(a;k!)\0) = (0\V r (a; k)V r (a; k')\0) = 



(A.4) 
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The operator product expansion of left handed vertex operators reads: 

a 2 

,V l{ -a- lk) V l{ a;k > = Y ^^) 
x (yV~ k\k\ a2 - 2 0(k)6(k' - k)a\ k ,_ k 

+Vk~¥\k'\ a2 - 2 6(k')6{k - fc'Hfc-fc, + ...) (A.5) 

*Vi(a;k)Vi(-a;k')* = - * Vft—a; -k)Vi(a; -k') * (A.6) 
Matrix elements of vertex operators between number states are: 

/ CY I1tt \ n '' k p- ikn l,kX 

<n|^(a;x)|0> = U k>0 -W^ (A.7) 



( a /2tt \" r ' fc e 



HF r(a ;x)|0}=n, >0 ^y T J — (A.8) 
Appendix B. Flow equation for hf\(B) 

The flow equations for the coefficient hf k (B) can be obtained from the commutator 
of the generator with the second part of the ansatz. Only the terms linear in a;^ are 
kept, and they come from the operator product expansion of two vertex operators: 



4(27r) 4 £^(p;i?)n[V;£) 



x (*VKa;p)VK-a; k - p') * (0\V r (-a;p)V r (a; -p')\0) 
+ * Vi(-a;p)Vi(a;p' - k) * (0\V r (a; -p')V r (-a;p)\0)) (B.l) 
By using the expressions given in Appendix A one finds the flow equation: 

S l )fT}\ ,/o_\4 / t \ 2a 2 



dh\:> k (B) _ 4(2tt) 4 /LY 

r(a 2 )r(a 2 -l) V27ry 



x 



£ t»(p; S)n« (p; B)p 2 " 2 - 2 + ^ u(p; Bjfi® (p; i?)p" 2 |p - fc^ 

(B.2) 

which is in second order in the coupling constant go . The integration of this equation 
turns out to be easy by using the flow equation ([32|) for the parameter Q k \(B). At 
leading order in the coupling constant the solution reads: 



*2<b) = i « 4 * ^ 



- 



x 



T(a 2 )r(a 2 - 1) \2tt J 
Y J (^)Hp-,B)p^- 2 + Y J {^f{p-,B)p a \p~kr 2 ^ | (B.3) 

.p>0 p>k 



We get (for ka < 1): 

m - 1 - 

r(a 2 )r(a 2 - 1) V2tt 



p>0 
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where we have defined il k l '{p;oo) = Q^(p;oo)/Vk. A similiar derivation can be done 
for the right-handed operators yielding h^\{B) ~ 1 — hf k . 



Appendix C. Sum rule 



The canonical commutation relations must be fulfilled during the entire flow: 

[a Uk (B),a\ k ,(B)] = ±-5 Kk , (C.l) 

From this requirement one can derive a consistency condition for the flowing operator. 
Evaluating this condition in the ground state gives: 



<5fc,fc' = hf k 5k,k' — hl k Sk,k' 



(2tt) 4 /' /. * 



T(a 2 ) 2 \2tt 



£ ng(p; B)p a \p + kf' 1 - «2(P; V 2 (P kf- 1 (C.2) 

\P>0 p>/c / 

Using some straightforward algebra one can verify that this is indeed fulfilled for the 
flow equations derived in this paper. 



References 



[1] Deutsch J 1991 Phys. Rev. A 43 2046 

[2] Srednicki M 1994 Phys. Rev. E 50 888 

[3] Rigol M, Dunjko V and Olshanii M 2008 Nature 452 854 

[4] Cramer M, Dawson C, Eisert J and Osborne T 2008 Phys. Rev. Lett. 100 030602 

[5] Reimann P 2008 Phys. Rev. Lett. 101 190403 

[6] Moeckel M and Kehrein S 2008 Phys. Rev. Lett. 100 175702 

[7] Moeckel M and Kehrein S 2009 Ann. Phys. 324 2146 

[8] Hackl A, Roosen D, Kehrein S and Hofstetter W 2009 Phys. Rev . Lett. 102 196601 

[9] Hackl A, Vojta M and Kehrein S 2009 Preprint arXiv:0908.3647, to appear in Phys. Rev. B 

[10] Eckstein M and Kollar M 2008 Phys. Rev. Lett. 10 120404 

[11] Eckstein M, Kollar M and Werner P 2009 Preprmf larXiv:0904.0976l 

[12] Manmana S R, Wessel S, Noack R M and Muramatsu A 2007 Phys. Rev. Lett. 98 210405 

[13] Kollath C, Laeuchli A and Altmann E 2007 Phys. Rev. Lett. 98 180601 

[14] Barthel T and Schollwoeck U 2008 Phys. Rev. Lett. 100 100601 

[15] Kinoshita T, Wenger T and Weiss D S 2006 Nature 440 900 

[16] Rigol M, Dunjko V, Yurovsky V and Olshanii M 2007 Phys. Rev. Lett. 98 050405 

[17] Calabrese P and Cardy J 2007 J. Stat. Mech. P06008 

[18] Cazalilla M A 2006 Phys. Rev. Lett. 97 156403 

[19] Grciner M, Mandel O, Haensch T and Bloch I 2002 Nature 419 51 

[20] Zamolodchikov A B 1977 Pisma. Zh. Eksp. Teor. Fiz. 25 194 

[21] Hackl A and Kehrein S 2008 Phys. Rev. B 78 092303 

[22] Hackl A and Kehrein S 2009 J. Phys.: Condens. Matter 21 015601 

[23] Wegner F 1994 Ann. Phys. (Leipzig) 3 77 

[24] Kehrein S 2006 The Flow Equation Approach to Many-Particle Systems (Springer) 

[25] Giamarchi T 2003 Quantum Physics in One Dimension (Oxford, Clarendon Press) 

[26] Iucci A and Cazalilla M A 2009 Preprewt larXiv:0903. 12051 

[27] Gogolin A O, Nersesyan A A and Tsvelik A M 1998 Bosonization and Strongly Correlated 

Systems (Cambridge, Cambridge University Press) 

[28] Coleman S 1975 Phys. Rev. D 11 2088 

[29] Kosterlitz J M and Thouless D J 1973 J. Phys. C: Solid State Phys 6 1181 

[30] Mandelstam S 1975 Phys. Rev. D 11 3026 

[31] von Delft J and Schoeller H 1998 Ann. Phys. 7 225 

[32] Kehrein S 1999 Phys. Rev. Lett. 83 4914 

[33] Kehrein S 2001 Nucl. Phys. B 592 512 



Sudden interaction quench in the quantum sine- Gordon model 



18 



[34] Glazek S D and Wilson K G 1993 Phys. Rev. D 48 5863 
[35] Glazek S D and Wilson K G 1994 Phys. Rev. D 49 4214 

[36] Fermi E, Pasta J and Ulam S 1965 Collected Papers of Enrico Fermi: 266. Studies of Non 

Linear Problems, Document LA-1940 (May 1955) (The University of Chicago Press) 
[37] Gallavotti G (Editor) 2007 The Fermi-Pasta- Ulam Problem - A Status Report (Springer) 



